GLOBAL STRONG SOLUTION TO THE DENSITY-DEPENDENT 
INCOMPRESSIBLE FLOW OF LIQUID CRYSTALS 



XIAOLI LI AND DEHUA WANG 

Abstract. The initial-boundary value problem for the density-dependent incompress- 
ible flow of liquid crystals is studied in a three-dimensional bounded smooth domain. 
For the initial density away from vacuum, the existence and uniqueness is established for 
both the local strong solution with large initial data and the global strong solution with 
small data. It is also proved that when the strong solution exists, a weak solution with 
the same data must be equal to the unique strong solution. 



1. Introduction 

Liquid crystals are substances that exhibit a phase of matter that has properties between 
those of a conventional liquid and those of a solid crystal. A liquid crystal may flow like a 
liquid, but its molecules may be oriented in a crystal-like way. There are many different 
types of liquid crystal phases, which can be distinguished based on their different optical 
properties. One of the most common liquid crystal phases is the nematic, where the 
molecules have no positional order, but they have long-range orientational order. The 
three-dimensional density-dependent incompressible flow of nematic liquid crystals can be 
governed by the following system of partial differential equations ([H [TT j |23 | [26] ) : 

d t p + V • (pu) = 0, (1.1a) 

d t (pu) + V • (pu u) + VP = pAu - AV • (Vd Vd) , (1.1b) 

dfd + u- Vd = 7 (Ad + |Vd| 2 d) , (1.1c) 

V-u = 0, (l.ld) 

where p denotes the density, u£l 3 the velocity d G § 2 (the unit sphere in R 3 ) the unit- 
vector field that represents the macroscopic molecular orientations, Pel the pressure 
(including both the hydrostatic part and the induced elastic part from the orientation 
field arising from the incompressibility V • u = 0) ; they all depend on the spatial variable 
x£l 3 and the time variable t > 0. The positive constants p, A, 7 stand for viscosity, the 
competition between kinetic energy and potential energy, and microscopic elastic relax- 
ation time or the Deborah number for the molecular orientation field, respectively. The 
term AV • (Vd Vd) in the stress tensor represents the anisotropic feature of the system. 
We set p = A = 7 = 1 since their exact values do not play any role in our analysis. The 
symbol Vd0Vd denotes a matrix whose (i, j)-th entry is d Xi d-d x d for 1 < i,j < 3, and it 
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is easy to see that Vd© Vd = (Vd) T Vd, where (Vd) T denotes the transpose of the 3x3 
matrix Vd. System (jl.ip is a simplified version, but still retains most of the interesting 
mathematical properties (without destroying the basic nonlinear structure) of the original 
Ericksen-Leslie model ([El fl^l fl6| [T7\ [23"1 125j ) for the hydrodynamics of nematic liquid 
crystals; see [U [Hll EH E21 E7] for more discussions on the relations of the two models. 
Both the Ericksen-Leslie system and the simplified one (jl.ip describe the macroscopic con- 
tinuum time evolution of liquid crystal materials under the influence of both the velocity 
and the orientation of crystals which can be derived from the averaging/coarse graining 
of the directions of rod-like liquid crystal molecules. In particular, there is a force term in 
the u-system (|l.lb|) depending on d; the left-hand side of the d-system (|l.lc|) stands for 
the kinematic transport by the flow field, while the right-hand side represents the internal 
relaxation due to the elastic energy. In many situations, the flow velocity field does disturb 
the alignment of molecules, and in turn a change in the alignment will induce velocity. 

We consider the initial-boundary value problem of system (jl.ip in a bounded smooth 
domain Qct 3 with the initial condition: 

(p, u, d) | t=0 = (po, u , d ), x G tt, (1.2) 

and the boundary condition: 

(u,$,d)|«i=(0,0), (L3 ) 

where p : Q -)• R + , u : ft -4 M 3 , d : ft ->• S 2 are given with compatibility, and v denotes 
the outer unit-normal vector field on <9ft. The boundary condition implies non-slip on the 
boundary and no contribution to the surface forces from the director field d. Roughly 
speaking, (|l.ip is a coupling between the incompressible inhomogeneous Navier-Stokes 
equations and the transported flow of harmonic maps. In the homogeneous case p = 1, 
(jl.ip becomes the hydrodynamic flow system of incompressible liquid crystals. In a series 
of papers [26]-[29], Lin and Liu addressed both the regularity and existence of global weak 
solutions to the Leslie system of variable length, i.e. when the Dirichlet energy 

- / |Vd| 2 dx, d : ft -4 § 2 , 
2 Jn 

is replaced by the Ginzburg-Landau energy 

(-|Vd| 2 + ^ ±J^_)dx (e > 0), d:ft^R 3 . 

2 Ae z 

In particular, for any fixed e > 0, they [26] proved the global existence of weak solu- 
tions with large initial data under the assumption that uo G L 2 (ft),do G -ff^ft) with 

3 

do IsnG H*(dQ) in the two-dimensional and three-dimensional cases. The existence and 
uniqueness of global classical solution was also obtained if uo G H 1 ^), do G H 2 (Q) when 
the fluid viscosity p is large enough. The partial regularity of the weak solution was inves- 
tigated in [27] (and also [5l 1171 132] ). similar to the classical theorem by Caffarelli, Kohn, 
and Nirenburgh [3j on the Navier-Stokes equations that asserts the one-dimensional para- 
bolic Hausdorff measure of the singular set of any suitable weak solution is zero. With the 
Ginzburg-Landau penalty function, the global strong and weak solutions to the compress- 
ible flow of liquid crystals were obtained in pi l3i| [35j 158] . See also P HH EU [32j E7] for 
some related discussions. For the incompressible version of system (jl.ip with constant den- 
sity, Lin-Lin- Wang |30j established the existence of global weak solutions that are smooth 
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away from at most finitely many singular times in any bounded smooth domain of M 2 , and 
we [Mj proved the global existence of strong solution in a bounded smooth domain of R 3 . 
For the compressible version of system (jl.lh . the one-dimensional classical solution was 
obtained in [91 [TO], [36] , and the blowup criteria of strong solutions were studied in [19j 120] • 
In this paper, we are interested in the existence and uniqueness of global strong solution 
(p,u,P,d) of (TTJ) in W 1 - r (n) x W 2 ^{9.f x W^(Q) x W 3 '«(fi) 3 with 3 < q < r < 
oo while assuming in addition that the initial density is bounded away from zero. By 
a strong solution, we mean a quadruplet (p, u, P, d) satisfying (jl.ip almost everywhere 
with the initial-boundary conditions (jl.2p - (jl.3p . Our strategy is to consider the following 
auxiliary problem: 

r d tP + v-Vp = o, 

pd t u- Au + VP = -pv-Vv- V- ((Vf) T Vf), 
' d 4 d-Ad = -v-Vf+|Vf| 2 f, 
,V-u = 0, 

for some given v 6 M 3 and f 6fi 3 . One of the motivations for such a strategy is that the 
continuity equation (ll.laj) is the transport equation of p, (ll.lbj) is the evolutionary density- 
dependent incompressible Navier-Stokes equation with the source term — AV • (Vd Vd), 
while (jl.lcp is the parabolic system in terms of d, therefore we can use a result of the 
transport equation (cf. Proposition 13. ip . the maximal regularities of the parabolic equa- 
tions (cf. Theorem 13. ip and density-dependent Stokes equations (cf. Theorem I3.2p . We 
first use an iteration method to establish the local existence and uniqueness of strong so- 
lution with general initial data. Then we prove the global existence by establishing some 
global estimates under the condition that the initial data are small in some sense. As 
system (|l.ip contains the Navier-Stokes equations as a subsystem, one cannot expect in 
general any better results than those for the Navier-Stokes equations. The uniqueness of 
global weak solution is always an open problem. We shall prove that when the strong 
solution exists, all the global weak solutions must be equal to the unique strong solution, 
which is called the weak-strong uniqueness. Similar results were obtained by Danchin [7] 
for the density-dependent incompressible viscous fluids in a bounded domain of M 2 with 
C 2+E boundary. We shall establish our results in the spirit of [7J, while developing new 
estimates for the crystal orientation field. Due to the particular structure of the equations 
for the velocity, especially the strongly nonlinear term (Vd) T Ad in the u-system, it will 
be necessary to obtain more regularities for the crystal orientation field. By developing 
more novel and subtle estimates, we will be able to finally establish the global existence of 
strong solution and weak-strong uniqueness for the initial-boundary value problem (jl.ip - 
(|1.3p of the density-dependent incompressible flow of liquid crystals. The results of this 
paper generalize our early results in |24| for the incompressible case with constant density. 
The analysis in this paper is much more difficult and complicated than that in [24] due to 
the appearance of non-constant density. 

The rest of the paper is organized as follows. In Section 2, we state our main results 
on local and global existence of strong solution, as well as the weak-strong uniqueness. In 
Section 3, we recall a standard result for the transport equation, the maximal regularities 
for the non-homogeneous non-stationary Stokes operator and the parabolic operator, and 
also some L°° estimates in the spatial variable. In Section 4, we give the proof of the local 
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existence. In Section 5, we prove the global existence. Finally in Section 6, we show the 
weak-strong uniqueness. 



2. Main Results 

In this section, we state our main results. If k > is an integer and p > 1, we denote 
by W k ' p the set of functions in L p (f2) whose derivatives of up to order k belong to L p (f2). 
For T > and a function space X, denote by L p (0,T;X) the set of Bochner measurable 
X-valued time dependent functions / such that t — > \\f\\x belongs to L p (0,T), and the 
corresponding Lebesgue norm is denoted by || • H^ppf)- We will consider the solutions in 
the functional spaces defined below. 

Definition 2.1. For T > and 1 < p, q, r < oo, we denote by Mj,' q ' r the set of quadruplets 
(p, u, P, d) such that 

u G C([0, T};D 1 Aq V' P ) n L p (0, T; W 2 ' q {VL) n Wq'^Q)), d t u G L p (0, T; L q (n)), V • u = 0; 

d G C([Q,T};B^ ) )nU 3 (0,T;W 3 ' q (n)), 8 t d G L p (0,T; L q (Q)); 

p G C([0,T]; W 1,r (p))\ P G LP(0,T;W 1 '«(n)) and / Pdx = 0. 

If r = oo, then p belongs to L°°(0, T; VF 1 ' oo (0))nC(0x [0, T]) instead of C([0, T]; W 1 ' 00 ^)). 
The corresponding norm is denoted by || • \\m p '' 



P,1,r . 

1 T 



We remark that the condition j n P (fx = in Definition 2.1 holds automatically if we 
replace P by 

p -w\L pdx 

__ i-ip 

in (jl.ip . Also, in the above definition, the space -D^ p stands for some fractional domain 

of the Stokes operator in L q (cf. Section 2.3 in [7]). Roughly, the vector-fields of D A p 

are vectors which have 2 — | derivatives in L 9 , are divergence- free, and vanish on <9fL The 

3(1— i) 

Besov space (for definition, see |2j) B QtP p can be regarded as the interpolation space 
between L q and W 3,g , that is, 

p <r 

3 ( 1 ^ 1 ) 1 1 1 

Moreover, we note that B q>p p VF 1 ' 9 if p > |. By the embedding W ,q L°° as g > 3, 

3(1--) 

one has -B^p p )■ L°°, which will be used repeatedly in this paper. 

The local existence will be shown by using an iterative method, and if the initial data 
are sufficiently small in some suitable function spaces, the solution is indeed global in time. 
More precisely, our existence result reads as follows. 

Theorem 2.1. Let be a bounded smooth domain o/M 3 . Assume that | < p < oo, 3 < 

1 l- 1 :? 3(1-^) 

q < r < oo, and po G W ' r (f2) with po > p for some p > 0, u G L) A J , d G B q ^ p p . 
Then, 
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(1) there exists Tq > such that system (|l.ip with the initial-boundary conditions 
(|1.2p - (jl.3p has a unique strong solution (p,u,P,d) E M^' r with < p < p on 
nx'[0,To]. 

(2) Moreover, there exist a constant v > depending on p, q, r, p, p, A, 7, £1, and a 
constant 5 > depending only on p, q, r, such that if 

v 



|U || ! 1 p + ||d || 3(1-1) < 



holds for the initial data, then the initial-boundary value problem (|1.1[> - (|1.3|) has a 
unique strong solution (p, u, P, d) E M^' q ' r for all T > 0. Furthermore, denoting 
by Ai £/ie /irst eigenvalue of the Dirichlet-Laplace operator in £1, for some constant 
C depending on p, A and 7, we have the following inequality for all t E R + : 

\\(^pu)(t)\\ L2 + ||Vd(i)|| z2 < Ce-^dl^uoll^ + HVdolb) (l + (^i*)a e ^ , 

with p = \\poWl 00 , and for some positive constant K depending only on \\pQ\\ w i, r , 
p, q, r, p, A, 7, p and tl, 

|| O, u, P, d) || M P,,,r<K || u 1| i.ip + Hdoll 3(i-i) • 



Similar to |30j . a weaA; solution to (jl.ip with the initial-boundary conditions (jl.2p - (|1.3p 
means a quadruplet (p, u, II, d) satisfying system (jl.ip in 17 x (0, T) for < T < 00 in the 
sense of distributions, i.e, 



/ / (pd t 4> + pu • V0) dxdt + / po0(-, 0) dx = 0, 
J Jnx(o,T) Jn 

/ / pu • (<9 t + u • X7(f>) dxdt + 11 Vu : V<p dx 

JJnx(o,T) JJnx(o,T) 

■ / p u ■ (/)(-, 0) dx + // (Vd Vd) : V(f> dxdt, 

Jn iisix(o,T) 



and 



d • d t 4> dxdt + 11 u • Vd • 4> dxdt +11 (Vd) : (V0) dxdt 

'Cx(0,T) J Jttx(0,T) JJnx(o,T) 

= / d o -0(-,O)dx + // |Vd| 2 d- 4> dxdt, 

Jn JJnx(o,T) 

for all <j) E C C °°(S1 x [0,T);R) or C c °°(0 x [0,T);IR 3 ). Moreover, (u,d) satisfies Cd in 
the sense of trace. In this weak formulation, the pressure II can be determined as in the 
Navier-Stokes equations (see [15]). 

Next, we will give a uniqueness result. Suppose 

peL°°(nx [0, T]) n C(0, T; U\VL)), V p > 1, 

a e l 2 '°°(o x [o,T])nw^°(n T ), pH 2 e r^Tji 1 ^)), 
d e L 00 ([o,T],H 1 (n))nL 2 ([o,T],H 2 (n)), vn e Lt(o,r ; Lt(fi)) 
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(for all T G (0, oo)) is a global weak solution to (|l.ip - (|1.3p . Then, we have the following 
energy inequality (cf. [30j Section 5 for the two-dimensional homogeneous case): 

~ I (p(i)|u(t)| 2 + |Vd(t)| 2 ) dx+ f [ (|Vu| 2 + |Ad+|Vd| 2 d| 2 ) dxdr 

* Jn Jo Jn in -<\ 

If ' 

< - / (/9 |u | 2 + |Vd | 2 ) dx, 
1 Jn 

for all t G (0,oo). 

However, as for the standard Navier-Stokes equations, the question of uniqueness in 
the above class of solutions remains open. For the same initial-boundary conditions, the 
relation between weak solutions and strong solutions can be formulated as: 

Theorem 2.2. Let Q,p,q,r be as in Theorem \2.1\ and /90) u (bdo satisfy the assumptions 
of Theorem \2.1\ Then any weak solution to (jl.ip - (jl,3p in the above class is unique and 
indeed is equal to its unique stronq solution. 

Usually, the uniqueness in the above theorem is called weak-strong uniqueness. For the 
similar results on the compressible Navier-Stokes equation, we refer the readers to |12|I31|. 



3. Maximal Regularity 

In this section, we recall a quite standard result for the transport equation and the 
maximal regularities for the parabolic operator and the non-homogeneous non-stationary 
Stokes operator, and prove some L°° estimates as well. 

For T > 0, 1 < p, q < oo, denote 

VM0,T) := (w^(o,T-,L^n))) 3 n(LP(o,T-,w^(n))) 3 . 

Throughout this paper, C stands for a generic positive constant. 

We first recall a result for the transport equation (cf. Proposition 3.1 in [7]): 

Proposition 3.1. Let Q be a Lipschitz domain of M 3 and v G (L 1 (0, T; Lip)) 3 be a 
solenoidal vector-field such that v • n = on dQ,. Let po £ W ' r (f2) with r G [1, oo]. Then 
the system 

(d t p + v-Vp = 0, 
\p\t=o = Po 

has a unique solution in L°°(0, T; W l,oc (fl,)) n C([0, T];r) q<00 W 1,q (£l)) if r = oo, and in 
C([0, T]; W l,r (£l)) if r < oo. Moreover, the following estimate holds: 

\\p(t)\\w^ < efo ^ T ^°° dT \\ Po \\ w i,r, t G [0,T]. 
If in addition p belongs to L p (Vt) for some p G [1, oo], then 

I lp = HajIIlp, t £ [0,T]. 



We first recall the maximal regularity for the parabolic operator (cf. Theorem 4.10.7 
and Remark 4.10.9 in [I]): 
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Theorem 3.1. Given 1 < p,q < oo, ujq G B^'^ and f G (L p (0, T; L 9 (1R 3 ))) 3 , the 
Cauchy problem 

(uJt - Aw = /, 
/ias a unique solution u) G Wg iP (0,T), and 



ll^llw,, p (0,T) < C I ||/||i^(i9) + |M| 3(1-1) j , 

where C is independent of ujq, f and T. Moreover, there exists a positive constant cq 
independent of f and T such that 



\M\W qp (0,T) > CO SUp \\U){t)\\ 3(1 _Xy 

te(o,T) B QtP p 

Now we recall the existence theorem (cf. Theorem 3.7 in [7]) for the linear system 
' pd t u- pAu + VP = f, J Q Pdx = 0, 

V-u = 0, (3.1) 
u| i=0 = u , u| 9n = 0. 

Theorem 3.2. Let O C R 3 a bounded domain with C 2+£ boundary, 1 < p < oo and 

3 < q < r < oo. Let u £ D A p ' and f G (i7(0, T; . ,4 ssume that the density p 
satisfies 

< p< p(x, t) < p < oo, (x, t) G O x (0, T), 
and /or some /3 G (0, 1], 

p G l°°(o, r ; w^cn)) n c^([o, r] ; L°°(n)). 

Then the system (|3.ip /ias a unique solution (u, P) suc/i //ia/ 

u g c([o,t\- d 1 ^' 13 ) n LP(o,T-w 2 ' q (n) nw^in)), 
s,u6j;(o,r;i«(fl)), 

and 

PGL p (0,T;^ ? (fi)). 

Moreover, there exists some constant C depending on p, g, r and VL such that for all t G 
[0, T], the following inequalities hold: 

i i_i 

pvp p||u(t)|| i-i, P + /*Hlz?(W«».«) + /5||<9 t u|| i P (iq) + \\P\\ L p {w i, q) 
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and 

ppi/~p\\u(t)\\^_x >p + \\pd t u,nV 2 u, VP\\ L v (Lq) 



< C ( e^+^)(^M 1- "l|u || + \\f\\ L > m ) + 

V D A q 



(3.3) 



where d(Q) is the diameter ofQ, £ p := /S//5, and 

B p (t) :=l + d(fi) (p^HVplUjo^r) 



2q 



M^(i) := ^Hplloj^nxip,*])' 
and i/ie exponents q,r* are numbers determined byp,q,r. 

Remark 3.1. The reader can also refer to Theorem 3.7 in [7J for more details about Theorem 



We notice that (|3.2p and (|3.3j) do not include the estimate for ||u||_£P ^ L qy Indeed, since 
we consider only in a bounded domain fi, then there exists a constant C = C(q,d(£l)) 
such that 

\\u\\ W 2, q = ||V 2 u|| L g + d(fi) _1 ||Vu|| L , + d(ny 2 \\u\\ Lq < C||V 2 u|| L9 , 

whenever u e PF 2 '«(fi) n W 1,9 (O) (cf. Proposition 2.4 in [7]). 

We also have the following two interpolation inequalities for the L°° estimates in the 
spatial variable (cf. Lemma 4.1 in [7J, also Lemmas 3.1,3.3 in [23]) which are useful in our 
proof. 

Lemma 3.1. Let 1 < p, q, r, s < oo satisfy 

p 3p 111 

2 2q s r q 

Then the following inequalities hold: 

IIV/IIj^oc) < CtMh/H 1 ^ ^ ||/||« (w2 , q) , 
liv/ll^^cri-^ii/n 1 ^ ||/||I, (W2>S)J 

for some constant C depending only on Q,p,q,r, and 

1-8 1 3 



p 2 2q 

Similarly, we can prove 

Lemma 3.2. Let 1 < p, q, r, s < oo satisfy 

2p p 111 

< -f - ^ < 1, - = - + -• 
6 q s r q 
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Then 



l|V/|| L , (i o 0) <C7T3-«||/|| 1 - e H/llWs,,), 



l|V/|L, (z .) < CT 



3 <j 



3( i_i) IU ll£P(W3.«)> 



for some constant C depending only on £l,p,q,r, and 

1-8 2 1 



P 3 q 

Proof. The proof is based on the applications of embedding and interpolation results in 
[2J. First, we notice that, from Theorem 6.4.5 in [2], 

3 _ 3_ 3 3 _3 1 — ft 2 1 
{Boo,So q , Boo,&)o,l = l With = - - 

p o q 

and the embedding (Theorem 6.2.4 in [2]): B^ 1 VF 1 ' 00 , we get 

||V/||zoo<q|/||%_3 11/11^3 3 . (3-4) 

R 9 n P q 

u oa,oo £>oo,oo 

We remark that 

B 3 ~,t«, W 3 > q ^ B* )00 B 3 ~,l (3.5) 

(cf. Theorems 6.5.1 and 6.2.4 in [2]). Therefore, according to (|3.4p . (|3.5p and by Holder's 
inequality, we deduce that 

IIV/II^) < C ( f \\f\fl:l\ WfK^dt) V < C ( f T \\f\\ P %%\\f\C, q dt) " 

<CTH||/f-0 11/11^^). 

rco/n 1 P ' \ J v y 

The proof of the second inequality is based on the fact that 

2-3-2 2 _3 1 — ft 2 1 

B r °! = (5 r / q ,B r>r q ) 0A ^L r with = o~- 

p 3 g 

(cf. Theorem 6.4.5 and Theorem 6.2.4 in [Sj), and that 

w 2 ' s # r 2 ; « , # s 2 ; > # 2 P ^ » , 

(cf. Theorems 6.2.4 and 6.5.1 in [2]). In fact, by Holder's inequality, we have 

iiv/ii^) < c ( [ T nv/f tf 3 nv/ir e 2 2 dt) < c ( / T ii v/ir^n v/i^ a ..dt ) P 

2 1 , „ 

< CT3"9 I 1 *" X ~ 6 



3(1-1) M -/ ll£f,(w 3 .«)- 

□ 
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Lemma 3.3. Let 1 < p, q, r, s < oo satisfy 

p p 111 
< | - - < 1, - = - + -• 
6 q s r q 



Then 



IIV 2 / || L p (LOO) < CT^ 



i i 



3 q 



3 {X _i) \\J IlLP (p,,), 



l|V 2 /|U {L ,) < CT 



3 q 



'TV-.- _ 3(1-1) ll-MIX*^ 3 ")' 



/or some constant C depending only on £l,p,q,r, and 

1-8 1 1 



p 3 q 

Proof. First, we notice that (cf. Theorem 6.4.5 in [2]) 

2_3_3 2—- 1 1 1 

{Boo,So 9 ,-Boo,i>)e,i = -8^,1 with = . 

p 6 q 

Hence, 

||v 2 /IIl» < c\\vf\\ w ^ < C||v/|| Bl < c||v/f 2 _3 liv/n 1 -^ 

oo.l D 9 p P 

^00,00 -000,00 

We remark that (cf. Theorems 6.2.4 and 6.5.1 in [2]) 

R P <^_v R P 9 W 2 ' 1 ? R 2 <^_v R 9 
9iP JJoo, oo j " * -"g.oo " u oo,oo- 

Thus, according to f|3.6[) and by applying Holder's inequality, we deduce that 

iiv 2 /ii l p (LOO) < c ( / T iiv/ir e 2 „| iiwirtf^ * ) P 

y ^00,00 $00,00 I 

<c( [ T \\f\\ p ^\\f\\ p %%dt] P 



3(1-1) \\J llLP(W3.9)- 

The proof of the second inequality is based on the fact that 

5^ = (5 r 1 ^~%s r 1 ;> 1 ^L r with — = 1-1 

p 3 q 

(cf. Theorem 6.4.5 and Theorem 6.2.4 in [2]), and that 
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(cf. Theorems 6.2.4 and 6.5.1 in [2]). In fact, by Holder's inequality, we have 

Iiv/IL^) < c ( f uv/ft? 3 \Wf\\ p \ sdt) p 

<c( [ T \wfi3\m\C*dt\ P 



L v T {w 2 ^y 



□ 



4. Local Existence 

In this section, we prove existence and uniqueness of strong solution on a short time 
interval, i.e., the local strong solution in Theorem 12.11 The proof will be divided into 
several steps, including constructing the approximate solutions by iteration, obtaining the 
uniform estimates, showing the convergence, consistency and uniqueness. 

4.1. Construction of approximate solutions. We initialize the construction of ap- 
proximate solutions by setting p° := po ; u ° := u o an d d° := do- Given (p n , u n , P n , d n ), 
Proposition 13. 1( Theorem 13.11 and Theorem 13.21 enable us to define respectively / o n+1 (x, t) 
as the (global) solution of the transport equation 



d t p n+1 + u n • Vp n+1 = 0, 



„n+l I „ 

P \t=o— po, 



(4.1) 



(4.2) 



d n+1 (x, t) as the (global) solution of 

jd t d n+1 - Ad n+1 = -u n • Vd n + |Vd n | 2 d n , 

jd^ 1 |t=o=d , a v d n + 1 \ da = o, 

and ((u n+1 (x,t),P n+1 (x,t)) as the (global) solution of 

' p n+1 d t u n+1 - Au n+1 + VP n+1 = -p n+1 u n ■ Vu n - V • ((Vd n ) T Vd n ), 
V- u n+1 = 0, J n P n+1 dx = 0, (4.3) 

„n+l| „ ,,n+l| n 

u T \ t =o — uo, u | a n — u. 

An argument by induction yields a sequence {(p n , u n , P n , d n )} ng N C Mj! q ' for all T > 0. 

4.2. Uniform estimate for some small fixed time T*. We aim at finding a positive 
time T* independent of n for which {(p n , u n , P n , d n )} ng N is uniformly bounded in the 
space Mj,' q ' r . Applying Proposition 13. 1 1 to (|4.ip . we get 

||p n+1 (t)llw^ < eti\\ v » n ^°° dT \\p \\ wl , r . (4.4) 

and 

min/9 n+1 (x, t) = p := minpo( x ) and max/9 n+1 (x, i) = p := maxpo(x). (4-5) 
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Since dtp n+1 = — u n • V p n+1 , then by Holder's inequality, we have 

V +1 £C(R + ;L s (fl)) 
with s = ^L- ( s = g if r = oo), and for t > 0, 

In order to apply Theorem 13.21 to (14.30 . we need to prove for some (3 G (0, 1], p n+1 G 
C /3 ([0,T];L oo (O)). Actually, noticing that by interpolation between L°°(0, T; W 1,r (fi)) 

and W 1 » oo (0,r;i»(n)), p n+l belongs to C^([0, T]; L°°(n)) whenever /3 G (0,^1) and it 

holds that 

||p n+1 || cf(LOO) < c(\\p n+1 \\ LT{m , r) + • (4.7) 

Here we have used Young's inequality. 

Hence, applying Theorem 13.21 to (|4.3p yields 

\\u n+1 (t)\\ 4 , + \\u n+1 \\ L P {w2 , q) + \\d t u n+1 \\ L P {Lq) + ||P n+1 || L?(iy i.) 
< Ce ct ^ ^||uo|| ol _i, p + ||u" • Vu" + V • ((Vd") T Vd")|| L?(i9) ^ , 



where 



m = (i + iip n+1 iu r( ^)) 70 (i + \\ P n+1 ^ 



for some positive exponent 70 depending only on p, q, r, /3, and the constant C depending 
only on p, q, r, p, p, Q, (3. Using (|4.6|) and (|4.7|) . we get 

m < C (l + ||p n+1 ||L r (H/^)) ,J (l + ||u n || L? >(W))* , 

where a depends only on p,q,r and (3. Therefore, 

||u" +1 (t)|| 4 , f + ||u" +1 || Lf(w , 2 , q) + ll^u^H^^) + \\P n+l \\ L p {w ,, q) 

D A q 

< ae C<i+l|p» +1 ll, r(wl>n ) CT (l+IKII^^))^|| Uo || i _ i + [| u » . Vu"|U M) (4.8) 



+ ||V • ((Vd n ) T Vd n> 
Applying Theorem 13.11 to (|4.2p . we obtain 

||d" +1 (t)|| 3(1 _ i) +||d^ 1 || W ^ (0i4) 



R 1 P' 
^ q ,p 



< 

Define 



C (ll d oll B 3 (1 -i) + II " u n • Vd™ + |Vd"| 2 d"|| L?(L9) ) . 



(4.9) 



U n (t) :=\\u n \\ + \\u n \\ L P {w2 , q) + \\d t u n \\ L P {Lq) 

L?{D Aq e' P ) 



+ \\d n \\ „ ri +||d n 



^d-i)^ ^ l|U IIW,. P (0,t). 
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U° = ||u || + ||d || 3(1-1), 

Q n {t) := \\p n \\ L f{w^) and g := \\p \\ w i, r . 
Hence, from (|4.4|) . (|4.8|) and (|4.9|) . we have the estimates for different cases as follows. 

Case 1. I — \ < i using Lemmas I3.1H3.31 we get, 

o q p 



|U 



<- r< C ' < ( 1+ " p ™ +1 "i?°(w 1 ' r )) ( 1 +!l un ll-c?°(i.9)) ? /||„ II I II., n || llY7,, n ll 

< ^ e V < >J V * 7 (ll u o|| i-l, p + l|U ||z, t °°(L9)l|VU ||l,P(Loo) 



+ l|Vd n ||_ L o ( iS )||Ad ri || i P( L oo)J, 
\d n+1 (t)\\ 3(1 _i ) + ||d" +1 || VV9!p(0it) 



B ' p 



< ^Ml do || 3(1 _i) + ||u n || L oo (L9) ||\7d n || i P (L oo) + ||d n || L c>o( L oo)||\7d n || L oo( i9) ||Vd n || L P( L oo)^, 
and 



+ tW(£p(i)) 3 ). 



(4.11) 



Ct^"p"^[/ n (t) < In 2, (4.12) 



Assuming that i is sufficiently small so that 
we get from (|4.1(jp that 

g n+1 (t) <2 6 o, (4.13) 

and from (|4.11[) that 

U n+1 (t) < Ce 2"C*(i+i/»(t)) *(i +eo r (U° + (t*-% + p-i + t^)(U n (t)) 2 

+t 2 3-i(u n {t)) 3 y 

9 1 1 3(1-1) 

Case S. | — g = p) then B q<p p W '°° with oo— denoting any positive number 
large enough (but not oo), using Lemmas 13.11 and 13.31 we get 

||d" +1 (t)|| 3(1-1) + l|d" +1 ||w„ p (0,*) 

Bq,p 

< C(l|do|| 3(1-1) + tP ||u n || i o ( i qr + )||Vd n || i oo( it x)--j 

+ * p l|d n |U°° (L°=)||Vd n || i oo( i , + )||Vd n || i oo( i oo-)^ , 
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and 

U n+l (t) < Ce ct ( 1+ e n+1 vY( 1+un Vy (U° + + t 1 - l «+t 1 ?){U n (t)Y 

+tp{u n (t)) a y 

Assuming that t is sufficiently small as that is in (|4.12p . we get f|4. 13j) and 

U n+l {t) 



< Ce 2 ' 



(4.14) 



< C(J|do|| 3{1 _i) + tp ||u n || L oo (L „)||Vd n || L c>o( i0 o) 

+ tp ll dn |lL t ° (L°°)ll Vdn ||L t °°(L9)ll Vdn ||Lf 1 (L °)))) 



3(1— i) 

Case 3. \ — \>h>\ — ^, then B q:P p )■ W 1 ' 00 , using Lemmas [3TT1 and EP1 we get 
||d" +1 (t)ll 3 (1 -i) + l|d n+1 Hw g , p (o,i) 

Bq,p 
3 



and again, by choosing t sufficiently small as that is in (|4.12j) . (|4.13j) and (I4.14D follow. 

Case 4. 1-^ = 1, then d]^J' P ^ W 1 ' 00 - and B^~^ ^ W 1 ' 00 , using Lemma E31 
we get, 

Q n+ \t)<e a e ctl '" u ^\ (4-15) 

and 

+ ||u n+1 ||z,p(h/2,9) + ||5ju n+ + ll- pn+1 |lLf(H/ 1 >'J) 



c^i+iip^n ^^i^y (i+\\u-w Lr{Lq) y / I 

S v t y >j \ t j M|Uq|| ^i^ + tfHu || Vu ||z,°°(z,°o-) 

+ ||Vd«|| ir(i?) ||Ad"|| i?(i o 0) ) ; 
U n+l {t) 

<Ce c < 1+ ^ +1 ^y( l+un ^ (u° + (tH +2ti-){U n (t)f + £{U n (t)f) . (4 ' 16) 
Assuming that i is sufficiently small so that 

Ct^p^it) <ln2, (4.17) 
we get from (|4.15p that (|4.13p holds and 

U n+1 (t) < Ce ^ct[^ n W)h^ S or(uO + (t H +2 tl)(C/"(t)) 2 + ^(L7™(t)) 3 ). (4.18) 
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Case 5. 1-±>1>1-1, then ^a/* ^ wl '°° and B w~^ ^ W l >°° , using 
Lemma [231 we get, 

||u" +1 (t)|| + ||u" +1 || L?(H/2 , 9) + \\d t vr +l \\ mLq) + \\p n+l \\ mw ^ q) 

D A q 

r ct(iM\ P ^ Lr l )\i+\\^\\ Lr{Lq) y / _i_rtiiii»n ii\7«»n 

< Oe v * ;/ v * ^ (Jl u o|| i_i, p + t p ||u ||x,oo ( W )||Vu ||i=°(i«>) 



+ ||Vdl L ~ (L , ) ||Adl L?(Loo) 

and (|4.16p follows. Moreover, by choosing t sufficiently small as in (|4.17p . (|4.13p and (14,181) 
follow. i 

Case 6. §-- = -, then D A P * ^ W 1 ' 00 and Bg^ p) ^ W 2 >°°- , we get 



u 



n+l '^ + ||u- +1 || if(w2 , s) + ||a t u" +1 || Lf(L9) + ||P" +1 || Lf(iyl , 9) 



Ci(l+||p n + 1 || i oo w l,r ) ) ( l+[|u n ||£,oo f Lq) ) 1 , n ,, , „, 

+ t^||Vd Jl || Lr{L9+) ||Ad™|| L P (ioo - 



and 

[T+H^Ce ^^ 1 '")'! 1 ^'^ (^° + 3^(^(t)) 2 + t^(C/"(t)) 3 ). (4-19) 
Assuming that i is sufficiently small as in (|4,17p . we get ()4.13j) and 

U n+ \t) < Ce 2 ° c < l+un{t ^ {l+e ^ (U° + 3i? (l^)) 2 + (C/"(t)) 3 ) . (4-20) 



Case 7. I - i > i then d\ P ' P ^ W 1 ' 00 and B?£ p) ^ W 2 '°°, we get 



\u n+1 (t)\\ i-I )P + l|u" +1 || L?(W / 2 , 9 ) + ||3 t u" +1 || L?(L9) + ||P n+1 || L? (^) 

i 



+ tp||Vd n || Lr(L9) ||Ad™|| Lf(Loo) j, 

and (|4.19p follows. Assuming that t is sufficiently small as in (|4.17p . then (|4.13p and (14.201) 
hold. 

Hence, for Cases 1, 2, and 3, if we assume that U n {t) < ACU° on [0,T*] with 



T* =min 



3pq — 2q — 3p 

/ In 2 \ 2 pi In 2 



U^V ' 2*C(1 + £ ) CT (1 + 4CC°)? ' 2 

i < 1 



1 \ 9-3 



16C 2 U°(3 + ACU°) 
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or 



1 < = min 



( 3pq-2q-3p 
In 2 \ ^ In 2 

4C2JJ ° J ' 2"C(1 + QoY (1 + 4CC/°) 



i > 
3 



max{p,^}' 



(4.22) 



WC 2 U°(3 + 4CU°), 
and, for Cases 4, 5, 6, and 7, if we assume that U n (t) < ACU° on [0,T*] with 



T* =min 



/ ln2 \FT In 2 

\iC 2 U°J Vc(l + o ) CT (l + 4CT/°p 



1 \max{p,^2_r 



(4.23) 



16C 2 U°(3 + 4CU°) 



< 1 



or 



1 < T* = min 



/ ln2 \FT In 2 

\AC 2 U°J Vc(l + o ) CT (1 + 4CT/°) 



j_ ■ 

/3 



1 



16C 2 [/°(3 + 4CTJ°) 



(4.24) 



then a direct computation yields 

U n+1 (t) < 4CU° on [0,7;]. 

Coming back to (|4.10p . we conclude that the sequence {(p n , u n , P n , d n )} is uniformly 
bounded in Mj,' q ' r . More precisely, we have proved the following estimates: 

Lemma 4.1. For all t £ [0,T*] with T* satisfying (j4.21|) or (|4.22p for Cases 1-3, and 
dP51) or (041 /or Cases .f 7, 

Q n (t) < 2g and U n {t) < ACU°. (4.25) 

4.3. Convergence of the approximate sequence. 

Lemma 4.2. There exists T$ such that {(p n , u n , P n , d n )},^ =1 is a Cauchy sequence in 
Mjf' r and thus converges. 

Proof. Let 

-n ._ p n+l _ p n^ Q n . = ^n+l _ u « pn ._ pn+1 _ pn ^ ._ jn+1 _ 

Define 

U n (t) :=||u"|| + ||u n || L P(vK2.») + \\dt^ n \\L p t {L=) 

+ l|VP n || L? ( L ,) + ||d n || 3(i-i) + ||d n || Wsp(0ii) . 

l T( b s,p ) 
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It is easy to verify that (p n , u n , P n , d n ) satisfies 
'd t p n + u n • Vp n = -u"' 1 • Vp n , 

p n+l d t u n - AH 11 + VP" = -p n (d t u n + u n • Vu") - p n (u n ■ Vu"^ 1 + u^ 1 • Vu r 

-V • ((Vd n_1 ) T Vd n ) - V • ((Vd^-^Vd"- 1 ) , 
d t d n - Ad n = -u"" 1 • Vd n - u"" 1 • Vd^ 1 + (Vd 71 ^™- 1 

+ ((Vd n + Vd™" 1 ) : Vd"- 1 ) d n ~ l , 
V-u n = 0, f n P n dx = 

with the initial-boundary conditions: 

(p n ,u n ,d n )\ t = = (0,0,0), 

(u B ,ft,d n )| fl o = (0 I 0). 

Applying Theorem 13.11 to 



n |2 An— 1 



<9 t d n - Ad n = -u"- 1 • Vd n - u"- 1 • Vd™- 1 + |Vd n |M 

+ ((Vd n + Vd™- 1 ) : Vd™- 1 ) d n ~\ 
^d"| i=0 = 0, ^d"| 9n = 0, 

and applying Theorem 13.21 to 

'p n+1 d t u n - Au" + VP n = -^{dtvC 1 + u n • Vu n ) - p n {u n ■ Vu"- 1 + u"- 1 • Vu n ~ 

-V • ((Vd n - x ) T Vd n ) - V • ((Vd n - 1 ) T Vd n - 1 ) , 

V • u n = 0, / n P n dx = 0, 
u n | t=0 = 0, u"| 9c = 0, 



we have 



U n (t) < C[\\p n {d t u n + u n • Vu n ) \\ l p {ls) + ||u n • Vu"" 1 + u"" 1 • Vu n ~ l \\ L 



+ 
+ 

+ 



< c(\\ff 



+ 
+ 
+ 
+ 
+ 

+ 



u 



n - 1 • Vd n + u n - 1 • Vd"" 1 ! 



v • ((vd"- 1 ) 1 vd") + v • ((Vd™- 1 ) 1 vd"- 1 ) \\ Ltm 

7~m-1 \ m — 1 ii 

Lf{L^){\\dtU n \\L p t {Li) + ll u "'ll-L| (-L«)l|Vu n || L P( Z/ oc,)) 
""■IIl^CL^IIVu™" 1 !!^^) + ||u n ~ 1 || i?C ( iS )||Vu n- ||if(ioo) 



| Vd n|2^n-l + ^ Vd n + Vd »-1) . yd"- 1 ) d"" 1 ) 



~- n — 111 i 1| 

\\V>{Lr) + II" I 

in— 1|| _ i 1\^-?ati— 1| 
in. — 1 II i 1 1 v7 j ri — 1| 



Vd n || L oo (Lq) ||Ad n \\iP(Lr) + l|Vd n 1 |U«>(i,s)||Ad n [| I/ P( L oo) 
Vd^lUco^llAd"- 1 ^^ + ||Vd-- 1 || ir(i?) ||Ad"- 1 || LP(i , ) 
u™ -1 II L t °°(L s ) II Vd^H^^oo) + ||u n " 1 || i oo( i , J )||Vd n ~ 1 || i P (L r) 
vdn 1 1 L c» (L , ) 1 1 Vd n 1 1 l p ( L r ) 1 1 d n ~ 1 1 1 L c» {LOO ) 

Vd^^lioo^^dlVd™!!^^) + ||Vd n_1 || L P (i0 o))||d n " 1 || L oo (L oo ) ^ 
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Note that if % — < | ; using Young's inequality, Lemma 13.11 yields 

II VU || L P (£r) < Ct* a« ||U || . +||U |U (W 2, S) ■ 

V ^°°(^a/ ) / 

1 '3 1 1 — iP 1 

If 2 — 2g > p> we have Z? As p W ' r so that the above inequality holds with the power 
of t replaced by £p. As for | — ^ = |, since 

II, ,n V7i^n— 1|| ^ ||,, n|| llY7,"^ n— 111 

ll u • Vu Wl p (l°) - H u \\l^(li+)\\\ /u \\l p (l^-) 

X— i p 1—-,P 

and D A v ' L q+ , D A p ' W ' r ~ with q+ (resp. r— ) slightly greater (resp. smaller) 
than q (resp. r), we still have 

\\u n ■Vu n - 1 \\ LP(LS) <tp\\u n \\ , i Jln^ii !_!„. 

The other terms such as HAd™ -1 ^^,.), || Vd n_1 H^p^r) and || Ad"!!^^) may be handled 

3(l-i) 

via the similar technique by using the Besov space B s , p p ■ 

To simplify the presentation, assume from now that | — | < | so that Lemmas I3.1H3.3I 

i 

can be applied. Otherwise, according to the arguments above, we would get tp instead of 

2_1 1__3_ 1_1 

i 3 i , t 2 2 i , or t 3 i below once their exponents are greater than or equal to -. 
Hence, for all t £ [0, T*], taking advantage of (|4.25p . the embedding 



3(i~i) ls qr 3 

B sp ' as s = — — > - 

q + r 2 



and Lemmas I3.HI3.31 we get 



c/ n (i) <c( ||p"|| Lr(i ,) + IIVu"- 1 !)^) +^-^||u"- 1 || Lr(LS) + HAd"- 1 !!^,) 



+ t3 «||vd n ^li-^ + tl ^iiu™- 1 !!^^) + ||vd n 1 \\ L v {L r ) 



2 _ 1 

L°°(L°°) I" << 3 9 ||vu llLf(L s ) 

13 11 2 1 



+ t3-^||d n - 1 || ft x> f r 00 -, +t3-g||Vd n - 1 



(4.26) 



Moreover, multiplying 

d t p n + u n • Vp n = -u 11 " 1 • Vp n 

by |p n | r_2 p n and integrating over Q, using V • u n = and the zero boundary condition, 
by Holder's inequality, we have 

1 d 



j7 \\P \\L r 

r at 



-- [ u"-v(|p n r') dx- [ IpT-Vu"" 1 • vp n dx 

r Jn Jn 



< IIl 1 - ll u ' v / 5 lli r 



By Holder's inequality, (14. 25ft and the embedding 

W 2 ' s (fl) L°°(ft) as s > |, 
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we eventually obtain 

\\p n {t)\\Lr < (\\v^\r)-V P n {r)\\ L rdr 
Jo 

< t 1_ p||u rt_1 H-" IIV^II (4.27) 



Inserting (pL"2Tj) into (14T261) . we get for t G [0, T*j, 

U n {t) < C(t l ~p +tU +ts~; + t§-|)i/«- i (t). 
If we choose To G (0, T*] such that 

C(T~* + T '~^ + r o *~« + Tf lq ) < i (4.28) 

then {(/3 n , u ra , P n , d n )} is a Cauchy sequence in M^f ,r and thus converges in 

We remark here that the time of existence To depends (continuously) on the norms of 
the data, on the bound for the density, on the domain and on the regularity parameters. 

4.4. The limit is a solution. Let (p, u, P, d) G Mjf' r be the limit of the sequence 
{(p n ,u n ,P n ,d n )}™ =1 in M%f' r . Passing to the limit in (g3} and (P5l) yields 

P<p{*,t)<p, (x, t) G n x [0, T ] and p £ L°°(0,To;W 1,r (Q)), 
u G L°°(0, T ; D~} ,P ) n 1/(0, T ; iy 2 ><^)), 9 t u G 1/(0, T ; L«(fi)), 

d G L^^TojsJ^JnLP^rosW^Cn)), 3 t d G #(0, 

P G L p (0,T ;W rl ' 9 (n)). 

We claim all those nonlinear terms in (|4.ip (|4.2p (j4.3|) converge to their corresponding 
terms in (jl.ip almost everywhere in x (0, To). Indeed, for a := ^r^(= 2 J+ r ) » 

Hu^-Vp^-u-VplUoo^) 

< llu" - u||^o( L .)||Vp n+1 |U«(ir) + ||u|| L oo (is) ||V/9 n+1 - VpHroo^r) 

< C (£>0||u n - u|| M P,o,r- + ||u|| L oo ( L s)\\p n+1 - p\\ M P>°>r 



as n — > oo, 



|| p n+l atU n+l_ Mu ||^ (LQ) 

< l|p n+1 ||Lf?(L-)II^U n+1 - <9 t u|| L P o(is) + - p\\L^ o (Lr)\\d t u\\ L P^ Ls) 

< C ^o||u n - u\\ M P£,r + UStull^^Hp^ 1 - p|| M *wJ 
— > as n — )• oo, 
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and 

ll _n+l_.n w-.n V7..II 



||p»+iu"-Vu"-pu-Vu|| L p 



< llp n+1 -pIIl-^II^II^^IIVu-II^^oc) 

+ ||p||L^(L«)||u n |U=?(L.)l|Vu n - Vu|| L P (Lr) 

+ l|p|U-(L-)l|u n " ullioo^llVuHi^^r) 

< C\\p n+l -p||L~(LT-)||u n || £ ,oc (L ,)||u n || L P ) ( W .2, g ) 

+ C||p||L«(Loo ) ||u n || L oo (£ , ) ||Vu n - Vu|| L P o(Lr) 

+ l|p||L-(L-)l|u n - u|U«(L-)ll Vu lli^ ) (i r ) 



< 



c((U°) 2 \\p n+1 - p\\ M , r + U°T ^ HplUoc (Loo) ||u« - u| 

+ r 5"^||p|| L?o(L oc ) ||u n -u|| M P^||u|| M ^) 

— > as n — > oo, 
due to u n — > u and p n+l — > p in Mj^' r as n — > oo. Hence, 

u n • V/9 n+1 -> u • Vp in (L p (0,T ;L a (n))) 3 ; 
p n+i dtU n+i pdtU in ( L P(o,r ;L Q (0))) 3 ; 
^+i u " . Vu n -> pu ■ Vu in (L p (0,T ;L Q (O))) 3 . 



Ml 



,s,r 



T 



Meanwhile, 



K-vd^-u-vdii^^ 

< ||d n - diUco^iivd"!!^^) + ||u|| L oc (LS) ||Vd" - vd|| L P o(io 

< C (u°\\d n - d|| M?r + ||u|| L co (L .)||d n - d|| M?r ) 
— ^ 0, as n — > oo, 



||V • ((Vd") T Vd") - V • ((Vd) T Vd)|| L , o(LQ) 

= ||iv ((Vd n + Vd) : (Vd n - Vd)) + (Vd n - Vd) T Ad n + (Vd) T (Ad" - Ad)\\ L P (La) 
< c(||Vd" - Vd|| L oo (LS) ||Ad" + Ad^^ + ||V<P + Vd|| L??o(LS) ||Ad" - Ad\\ L P To{Lr) 
+ \\Vd n - Vd|| L oo (is) ||Ad"|| L , o(Lr) + ||Vd|| L??)(La) ||Ad™ - Ad\\ LPTo(Lr) ) 

<CT ( f""(C/ + ||d|| Mr „)||d--d|| M ^ 
— > 0, as n — > oo. 
Then, we have 

Vu n -Vd n ^u-Vd in (L p (0, T ; L s {VL))f- 
V • ((Vd n ) T Vd n ) -> V • ((Vd) T Vd) in (L p {0, T ; L a (n))f. 
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If s > 3, we have 



|||vd"| 2 d"-|vd| 2 d|| iPTo(L§) 

< uvdi^) (livd-n^^^iid" - d\\ L%{Li) + ildi^^iivd" - vd|| i§?o(L5) 

+ ||d|| L oc ( ^ ) ||Vd|| iPTo(i9) ||Vd" - Vd||^ (L§) 

< CHVdlicc^,) (\\Vd n \\ L ^ {w i, q) \\d n - d|Uoo (L .) + ||d|| L , o(wl , s)) ||Vd" - Vd|| L oc (LS) ) 
+ C||d|| £ oo (W ri,. )) ||Vd|| £?b(£ r ) ||Vd n - Vd|| L oo (£ .) 

< CU° (U° + ||d|| M|r ) ||d" - d|| M?;r + CT* "*||d||^ r ||d" - d\\ Kr 



0, as n — > oo. 



If | < s < 3, we have 



HI Vd "l 2dn -|vd| 2 d||^ (L|) 

< C (|||Vdf (d" - d)\\ L%{Lh + ||Vd» : (Vd" - Vd)d|| £J) 
+ ||Vd:(Vd»-Vd)d|I^ (L1&) 



< CUVdl^) (||Vd"|| £ , o(£ ^||d" - d||^ (L§) + Hdll^^cojIlVd" - Vd||^ (L§) ) 

+ l|d|l L ??o (L^) l|Vd|l ^o^) l|Vdn - Vd|l 4^) 

< C||Vd w || £ oc (£9) (llVd^H^^Hd" - d|| £ ~ (£s) + ||d|| i? , (w2 , 3)) ||Vd n - Vd|| £ oc (£s) ) 
+ C||d|| £ oc (w i, s)) ||Vd|| L P o(w2 , s) ||Vd" - Vd|| £ =o (£ .) 

< CU° (U° + Hdll^-r) ||d" - dH^-r + C||d|| 2 ,r\\d n ~ d\\^,r 
\ T J T Tn T 



TO 

0, as n — > oo. 
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The case s = 3 may be handled by noticing that we also have 

|||vd«| 2 d™-|vd| 2 d|| iPTo(L§) 

< livdi^ (||Vd»[|^)||d» - d\\ L%[Lh + Hdii^^iivd- - vd|| L??o{ii) ) 

+ ||d|| ^^iVdH^^jHVd" - VdH^^j 

< C||Vd n || L? o( i9) (j|Vd n || L ^ (H/M) ||d n - d|| L oo (L 3) + ||d|| L P (W 2,3)) || Vd n - Vd|| L oo (L 3) 
+ ^l|d||L § ? ) ( H /i>3 ) ||Vd|| L? , o(w/2 ,3 ) ||Vd' 1 - Vd^oo^) 

< CU° (u° + ||d|| M ,3,) ||d» - d|| M? ,3„ + C||d||^ r ||d» - dWug.r 

— > 0, as n — > oo. 
Hence, we finally get 

|Vd n | 2 d n -> |Vd| 2 d in (L P (0,T ;L§(^))) 3 . 

Thus, passing to the limit in (jITTJ) . (f4T2|) and P~3|) as n -> oo , since L s (0) ^ L Q (f2) M> 

Lf(O), we conclude that (11. ip holds in (L p (0,To; L§ (O))) 3 and therefore almost every- 
where in f2 x (0, To). 

Multiply the d-system (jl.lcp by d, we obtain 

^(|d| 2 ) + i u • V(|d| 2 ) = Ad • d + |Vd| 2 |d| 2 . 

Since 

A(|d| 2 ) = 2|Vd| 2 + 2d-(Ad), 

then it follows that 

^(|d| 2 ) + l -u • V(|d| 2 ) = \ A(|d| 2 ) - |Vd| 2 + |Vd| 2 |d| 2 . 

Therefore, it is easy to deduce that 

9 t (|d| 2 - 1) - A(|d| 2 - 1) + u • V(|d| 2 - 1) - 2|Vd| 2 (|d| 2 - 1) = 0. (4.29) 

Multiplying (|4.29p by (|d| 2 — 1) and then integrating over Q, using (jl.ldp and (jl.3p . we 
get the following inequality: 



- f (|d| 2 - l) 2 dx < 4 I |Vd| 2 (|d| 2 - l) 2 dx 
dt Jn Jn 

<4||Vd||i« / (|d| 2 - l) 2 dx. 



(4.30) 



Remark that interpolation between L°°(0,T ; W l ' q (n)) and L p (0,T ; W 3 ' q (Q)) shows that 
for some positive a > ±, d belongs to L 2 (0, T ; H 2+a (Q)) and that llVd]] 2 ^ £ L^O.Tb). 
Notice that 

/ (|d| 2 - l) 2 dx = 0, at time t = 0. 
Thus, using (|4.30p together with Gronwall's inequality, it yields |d| = 1 in O x (0, T ). 
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4.5. Uniqueness and continuity. Let (pi, Ui, Pi, di) and (p 2 , U2, P2, d 2 ) be two solu- 
tions to (jl.ip with the initial-boundary conditions (j 1 . 2 [) (|l,3p . Denote 

p = pi~P2, u = ui-u 2 , P = Pi-P 2 , d = di - d 2 . 
Note that the quadruplet (p, u, P, d) satisfies the following system: 

d t p + ui ■ Vp = -u • Vp 2 , 

pid t u - Au + VP = -p(d t u 2 + ui • Vui) - p 2 (u ■ Vui + u 2 • Vu) 

-V- ((Vdi) T Vd) - V- ((Vd) T Vd 2 ), 
d t d - Ad = -ui • Vd - u • Vd 2 + |Vdi| 2 d + ((Vdj + Vd 2 ) : Vd) d 2 , 
V- u = 0, j Q P dx = 

with the initial-boundary conditions: 

(p,u,d)| i=0 = (0,0,0), (u,d v d)\ dQ = (0,0). 

Using the same argument for p n in Subsection 4.3, for all t S [0,To], we have 

\\p{t)\W < /Vp2(T)|M|u(T)|| Lo =dT 







<t l " 11 V^||x,oo( Zr r) ||u||x^(x,«») ^ 4 ' 31 ^ 
1 1. 

< Ct p ||/52||Lf ) (W 1 >'-)ll u llLf , (W 2 > s )- 

On the one hand, since pi,p 2 G L°°(0, T ; W 1 ^^)) n W 1 ' 00 ^, T ; L S (Q)) implies that 

PI1P2 G C^([0, 7b]; L°°(f2)) whenever /3 6 (0, j—j), then Theorem 13.21 yields, for some 

i 

constant C depending on Tq,p, q, r, p, p, O, f3 and on the norm of p\ in L°°(0, To; W 1,r (Q))r\ 
C /3 (0,r ;L oo (fi)), and for all t G [0,T ], 

||u(i)|| + ||u||xP(w 2 .«) + \\dtu\\ L P( L s) + ||P||lP(^i,3) 

< C(\\p{d t \i2 + ui • Vui)|| L P (La) + ||p 2 u • Vui|| l p (Ls) + ||p 2 u 2 • Vu|| l p (Ls) 

+ || V • ((Vdx) T Vd) || Z?(LS) + || V • ((Vd) T Vd 2 ) \\ l p {ls) ) 

< C {\\p\\ L c a ( yL r- ) {\\dt\\2\\ L p t {L'i) + HuilUfCLajllVuiH^P^oo)) + ||Vui|| l p (loo) ||u|| L oo (LS ) 

+ ll u 2|kf (£9)||Vu|| L p (ir .) + ||Vdi|| L oo (L9) ||Ad|| L P (Lr:) + IIAdill^P^^lVdll^oo^) 
+ l|Vd 2 || L oo (L?) ||Ad|| L P (ir) + || Ad 2 ||£j(i«») llVdllxoo^ 



< C\\\p\\ L ^( yL r){\\d t \l2\\L P t {L'l) + ll U l|l 1-1, H U l|lLP(iy2>'J)) + l|ui|| L P(Vl/2,9)||u|| L oo( iS ) 

\ L^{D Aq ' J ) 

+ ||u 2 || !_i || Vu \\mLr) + (II di|| 3(1 _i) +||d 2 || 8(1 _i) )||Ad|U (ir) 



(4.32) 



+ (ll d l llLP(W3.9) + ll d 2||Lf'( H /3,, ) )||Vd|| L oo {La) j. 
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On the other hand, Theorem 13.11 yields, for some constant C independent of To, 
\\d(t)\\ 3(1 _£) + ||d||w s , p (o,i) 

< C (\\ - ui • Vd- u ■ Vd 2 + |Vdi| 2 d + ((Vdi + Vd 2 ) : Vd)d 2 || i p (L , ) 

< ^([lUlllif (L^IIVdll^P^r) + ||Vd2|| L P( L oo)||u|| L c>c (Ls) 

+ l|Vdl|| L oo( W ) || Vdl ||iP(ioo) ||d||£,oc( L r) 

+ (llVdiH^^oo) + ||Vd 2 || L?(Loo) )||Vd|| Lr(L , ) ||d 2 || Lr{L ^ ) ) ( 4 - 33 ) 

- c '(ll Ul ll r oo/n 1 -p^" va " £ ?^ p ) + ll d2 Hif(w a .«)ll a IUr(i') 

+ ll d llL w , B S(l-i),ll d l|lif(W».«)l|d|U|0(Lr) 

+ ll d 2|| ^ 3(i-i) (||di|| Z f(pF 2 .8) + ll d 2|lLf(iy 2 >9))||Vd|| i oo (is) 

We remark here that Holder's inequality and the embedding W l,q {Q) > (q > 3) 

have been employed repeatedly in both (|4.32p and (|4.33j) . 
Lemmas I3.1H3.3I yield, by use of Young's inequality, 

||Vu|| L P (Lr) < Ct5-*(||u|| i + ||u|| if(w2 , s) ), 

L t \ U A S ) 

||Vd|| if(ir) < Cti-i([|d[| aa-j). + PIIl?^)), 

L-t ( B s,p I 



\Ad\\ mLT) < Ct* *(||d|| 3(1-1) + ||d|| if(w3 , s) ). 



Define 



x (t) HMIjW) + || u || ^ x i p + ||u|| L P (H/2 , a) + ||0tu|| L P (LS) 

L t°( D A s ) 

+ l|VP|| Lf(Ls) + ||d|| 3(i-i) + ||d|| Wsij , (0 , t ). 

3(1 — l) 

Thus, combining (ITOjl - flPgD and B s , p p ^ W 1,s (fl) ^ L r (tt), we have 



X(t)<CU p\\P2\\l™(W^)( 1 + W d t u 2\\ L P(Li) + ll u l|l ^ i-i P ll u l \\l?(W*«)) 

L T\ D A q ) 

1___3_ 2_1 
+ £ 2 2 9||U 2 || !_1 +t 3 "||ui|| L _l + ||ui|| L P (W 2, 9) 

i t °°P A / ) iS°(^, p ) ' 

1_1 . . 
+ * 3 9 (ll d l|l 3(1-1) + ||d 2 || 3 (1 _ij )+ ^2 X1(^2,,) 

+ ll d i|lL£(iy 3 ><0 + ||d 2 || z ,p (Ty 3, 5 ) + ||di|| 3(1 _ij ||di|| L P (H/ 2, 9 ) 

L t°( B q,p ) 

+ ||d 2 || ^ 3(1 _i, (||di|| L P (W 2,,) + ||d 2 || L P (VK 2, 9) )|x(t). 
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Now, choosing r\ so small that the term between brackets is less than | for t = r\ enables 
us to conclude that X = on [0,??]. As the constant C does not depend on 77, a standard 
induction argument yields the uniqueness on [0, To]. 

Finally, as p satisfies a transport equation with data in W l,r (£l), u satisfies 

^u-Au + VPG (L p (0,t;L"(O))) 3 , 

and d satisfies 

d t d- Ad G (L p (0,t;L 9 (fi))) 3 , 
then, Proposition EU Theorems UJ insure that p £ C([0, T ]; W^O)) (if r / 00), 
u G C([0,T ];^/' P ) and d G C([0, T ]; B^"* 5 ). 

Remark 4.1. Following the argument of uniqueness and continuity, we can also easily prove 
that if (/>!, Ui, Pi, di) and (p2, u 2 , P2> ^2) are solutions to f|l.lj) - ([1.3p with different initial 
data (po) u o> dg) an< ^ (Po^O'^o)' then the following estimate holds true on [0, To]: 

\\p(t)\\ L r + ||u(t)|| + \\u\\ L P {W 2, a) + \\d t u\\ L P {LS) 

D A S 

+ ll-Pllifcw 1 .*) + P(*)ll 3(i-i) + ||d||w a ,p(o,t) 
< C llpolU*- + ||uo[| i-i.p + ||do|| 3(i-i) , 

\ "s,P / 

where />o := Pq ~ Pq-> uo := Uq — Uq, do := dj — dp. Combining with Theorem 12.11 
we conclude that for small enough T, the map (po,Uo,do) — > (p, u, P, d) is Lipschitz 

1-A,p 3(1-1) 

continuous from bounded sets of W ,r x Z)^ p x p to 

C([0, T];I/(fi)) x (C([0, T]; n (^(0, T; L s {Q))f n (L"(0, T; ^ 2 ' s (fi))) 3 ) 

x IJP^T-W^in)) x (c([0,T];P s 3 S" p) )n (W 1 *(0,T;L-(n))) 3 n (i7(0,T;T^ 2 ' s (O))) 3 ). 



5. Global Existence 

In this section, we prove that, if the initial data of velocity and orientation field is 
sufficiently small in appropriate norms, the local strong solution (p, u, P, d) of (|l.ip - (|1.3|) 
established in the previous section is indeed global in time. 

5.1. Estimates for ||u|| L 2 and ||d|| L 2. 

Lemma 5.1. Let fl,p,q,r be as in Theorem \2. 1\ and let (p, u,P, d) G M^ q ' T be a solution 
to (|l.ip - (jl.3p on £1 x [0,To]. Then the following inequality holds true for all t G [0, To] : 

||(^u)(i)||| 2 + ||Vd(i)|| 

< e"^ 1 * (IIV^u ||| 2 + ||Vd ||| 2 ) (l + ^j-te^ , 
where X\ stands for the first eigenvalue of the Dirichlet- Laplace operator in f2. 
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Proof. Due to the inhomogeneous incompressible character the flows we are dealing with, 
the natural framework in which we shall work is that of the solenoidal vector field of 
L 2 (S1) 3 . Note that 

u G C([0, T ];D^**) n (1/(0, T ; W 2 «(V) n * , 

d e cc^ro];^ 1 "*^ (LP(o,r o; w 3 '9(n))) 3 . 

And, since 

l- 1 :? 2(l-i) 

(see Proposition 2.5 in [7]) where 

X q = {z£ L q (nf | V • z = in n and z • n = on dtl}, 
then, when 1 < p < 2, by the standard interpolation inequality 

L°°(0, T ; L«(fi)) n 1/(0, T ; W 2 >«(n)) C L 2 (0, T ; # 1+a (ft)), 



where 



we have 



2 oo p p 2 3 

,3 



u € C([0, T ]; n (L 2 (0, T ; tf 1 ^)))' , (5.1) 

d G C([0,T ];# 1+Q (ft))n (L 2 (0,T o;J F/ 2+a (ft))) 3 . (5.2) 
When 2 < p < oo, USHJ-Q hold obviously due to T^ 2 '*?^) <^-> H 2 (Q) as g > 3. 

Now, p is continuous in (t,x), u G C([0, T ]; tf a (ft)) n (L 2 (0, T ; if 1+a (ft))) 3 and 

d G C{[0,T ]-H 1+a (n)) n (L 2 (0,T ;# 2+a (ft))) 3 . This enables us to justify the follow- 
ing computations. 

Taking the L 2 scalar product in (jl.lbl) with u and performing integration by parts, 
using the continuity equation (II. lap , we obtain 

55/ n "' u l 2dx + / n ' Vu ' 2dx 

= -- / |u| 2 V • (pu) dx - / pu ■ Vu • u dx - / u • ((Vd) T Ad) dx (5.3) 
2 Jn Jn Jn 



[ u ((Vd) T Ad) dx. 
Jn 



Here we have used the facts 

V • (Vd Vd) = V P^j^) + (Vd) T Ad, 
and V • u = in £1, u = on dQ, as well as 

J VP • u dx = J V • u dx = 0. 

Multiplying (|l.lb|) by —(Ad + |Vd| 2 d) and integrating over f2, we obtain 

- / <9 t d-Addx- / (u-Vd)-Addx = - f | Ad + |Vd| 2 d| 2 dx. 
Jn Jn Jn 
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Here we have used the fact that |d| = 1 to get 

{d t d + u • Vd) • |Vd| 2 d = - (|Vd| 2 di(|d| 2 ) + u • V|d| 2 |Vd| 2 ) = 0. 

Since d u d = on Oft, integrating by parts, we have 

/ d t d- Ad dx = -~ I |Vd| 2 dx. 
Jn 2dtJ n l 1 

Hence we obtain 

- — / |Vd| 2 dx + / |Ad + |Vd| 2 d| 2 dx = I (u • Vd) • Ad dx. (5.4) 
2 dt Jn Jn Jn 

By adding (|5.3|) and (|5.4|) . we eventually get the identity: 

\— I (p|u| 2 + |Vd| 2 ) dx+ [ (|Vu| 2 + | Ad + |Vd| 2 d| 2 ) dx = 0. (5.5) 
2 dt Jn Jn 

Since Vd G L 2 (0, T ; H 1 ^)) and |d| = 1, we have 

Ad-d + |Vd| 2 = 0, (5.6) 

and then 

f |Ad + |Vd| 2 d| 2 ) dx = [ |Adxd| 2 dx. 
Jn Jn 

Now, by virtue of the Poincare inequality ||Vu|| 2 2 > Ai||u|| 2 2 , we get 

\j t (llv^lli- + ll vd ll!0 + j\W~p*\\h < o> 



i.e. 



d ( Hi t|l _ ll2 \ , 2*L t d 



- ^^HVpuII^J +e ^ t -||Vd|| 2 :2 < 0. (5.7) 
Integrating (|5.7|) from to t, we obtain 

{\\(^pu)(t)\\h + m(t)\\h) 

2\, ft 2\ 1 

< ||V^u || 2 2 + HVdollia + - 1 / e p T ||Vd(r)|| 2 2 dr. 

P Jo 

It follows from Gronwall's inequality that 

2A 1+ 2A1 a u 

e^llVdH 2 , < (||^u || 2 2 + ||Vd || 2 2 )(l + -^te~\ 

and furthermore, 

e^dlVpuH^ + ||Vd||| 2 ) < (Hv^uollia + ||Vd ||£ 2 )(2 + ^ie^ - e> + ^i). 

□ 

Usually, (|5.5p is called the basic energy law governing the system (|l.lj) - (|1.3[) . It reflects 
the energy dissipation property of the flow of liquid crystals. 
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5.2. A more explicit lower bound for the existence time. We denote by T* the 
maximal existence time for (p, u, P, d) which means (p, u, P, d) cannot be continued be- 
yond T* into a strong solution of (jl.lj) - (|1.3|l . Let us first state a continuation criterion: 

Lemma 5.2. Let po,uo,do be as in Theorem \2.1\ and assume that system (jl.ip with the 
initial-boundary conditions f)1.2[) - f|1.3|) has a strong solution on a finite time interval [0, T*) 
with 

p€L oo (0,r*;W 1 ' r (fi)), inf p(x,t) > 0, 

t<T*, xen 

u € L oo (0,T*;D 1 A J' P ) and d G L°°(0, T*; P^p"^). 
T/ien (p, u, P, d) can 6e continued beyond T* into a strong solution of (|l.ip - (|1.3p . 

Proof. Indeed, a positive lower bound T for the existence time has already been obtained 
in the proof of Theorem 12.11 (see (|4.2ip(|4.22p and (|4.28p ) when (po,uo,do) remains in a 
bounded set of 

-, l--,P 3(l-i) 

W 1 ' 1 x D Ag p x B g " p p> 

with in addition inf xG Qpo(x) — P f° r a fixed p > 0. Hence system with initial 

density p(T* — initial velocity u(T* — ?) and initial orientation field d(T* — ^) has 
a unique strong solution on [0,T] which provides a continuation of the strong solution 
beyond T*. □ 

Combining Lemma 15. II and Lemma 15.21 enables us to get the following result: 

Proposition 5.1. Let po, uq, do be as in Theorem \2.1\ and let (p, u, P, d) denote the corre- 
sponding strong solution of fll.lji — fll.3j) . Then there exists some constant C depending on 
p, q, r, p, A, 7, £1 and p, such that, the maximal existence time T* for (p, u, P, d) satisfies 

C 

T ~ (I + HpoII^m-W ) 1 
for some positive exponents k and i depending only on the regularity parameters. 

Proof. Fix a T < T*. We aim at proving that if f < C(l + \\po\\w 1 ^)~ K ( u °)~'' for a 
convenient choice of C, k and i then (p, u, P, d) may be bounded in M^ q ' r by a function 
depending only on the data. Then Lemma 15.21 will entail Proposition 15.11 
Define 

+ l|VP|| L p (L9) + ||d|| 3(1 -l) + l|d||w 9 , P (o,t)- 
According to Theorems I3.HI3.2I and (|4.7p , we have 

G(t) < Cfe(t)(||u || + ||u • Vu|| Lf(i9) + ||V • ((Vd) T Vd)|| Lf(i9) ) 

+ C p (0||u|| l p (L9) + ||d || 3(i-i) + || - u • Vd + |Vd| 2 d|| L P (Lq 



(5.8) 



J 1-V 



where C = C{p, q, r, Q,, p, p). 
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Combining the Gagliardo-Nirenberg-Sobolev inequality and Young's inequality yields, 
for all e > 0, 

||u||l, < C(e\\u\\ W 2, q +e 1 -ff||u|| L2 ) with 6 = — ^— . (5.9) 

7q — 6 



We note that Lemma 15 . 1 1 insures 

HU r(L2) <C(l + th)U°, 
then employing Holder's inequality, we have 

H\L*(LZ)<th\"\\L~(L2)<Ctr(l + £)U . (5.10) 

Moreover, We get 

||u • Vu|| L?(Lg) < llulli-^HVull^^.) < Ct^G 2 (t), (5.11) 

||V • ((Vd) T Vd)|| if(i9) < CIlVdlUac^HAdlliP^oc) < Ct^G\t), (5.12) 

|| -u- Vd+ |Vd| 2 d|| L P (L9) < ||u||z|=(i9)||Vd||^(£oo) + ||Vd||i,«» (i ,)||Vd||£y(£oo) 

i 3 (5.13) 
< Ct*~^G 2 {t). 

Here we have used the fact that |d| = 1. Hence plugging (15.9p - (|5.13p in (15.81) while taking 
e = eC~ 1 (t) with e suitably small, we get 

G(t) < c(B 2 p {t){U° + {t^~^ +t*~«)G 2 (t)) +cl(t)tp{l + t^)U° + U° + t^'^G 2 (t)Y 

(5.14) 

On the other hand, using the same argument as for p k+1 in Subsection 4.2, we obtain 

UVplUoo^) < \\p \\ w i,re cti ~^ G ^, (5.15) 
\\p\\ C e(Lo°) ^ C {\\p\\l?(W^) + W^PhriL 3 )) 

< C(\\p\\ L oo iW l,r) + \\u\\ L oo( Lq) \\Vp\\ L oa( L r}) (5.16) 

<C||Vp|| Lr(L ,)(l + G(t)). 

Then, according to the definitions of B p {t) and C p (t) in Theorem 3.7 in [7J, using (|5.15p 
and (|5.16p . we eventually get, 

B p (t) < C7e ct§ "^* G W(l + \\p \\w^)^, (5.17) 

< Ce c ' ti "-"*^)((l + || >O0 || v ^ 1 ,.)' yi + (l + ||^ ||^ 1 -)' y2 |l/'o||^ a ^(l + G'(i))^), (5.18) 

where 71 and 72 depend only on p,q,r and (3. 

Plugging (15.17j) - (15.18p in (I5.14p . for some positive exponents 5\ and 62, we have 

G(t) <Ce cti ^ G ^(l + \\p \\ w i $r )^ (c/°(l + ti(l + ^)(l + G(t)) 52 ) 

+ (t*~^ +^ _ 9)G 2 (t)^). 

Assume that T has been chosen such that 

G(f)<8C(l + \\po\\ wl ,r) 5i U°. (5.19) 
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This is possible because of the continuity of the function 1 1— > G(t). Noticing that G(t) is 
increasing in t, then a standard induction argument shows (|5.19p is satisfied at time t < T 
with a strict inequality whenever the following three inequalities are satisfied: 

8C 2 (1 + \\po\\w^) 5l U°t^'p~^ < ln2, 

(l + 8C(l + \\po\\ w i,r)^U ) S2 tl(l + t^ < 1, 

64C 2 (1 + \\po\\ w i,r) 2Sl U°(^'^ + iH) < 2. 

Hence Lemma 15.21 enables us to continue the solution (p, u, P, d) beyond T. 

The proof of Proposition 15.11 is complete. □ 

5.3. The case of a small initial velocity and orientation field. Proposition 15.11 
insures that the existence time of a strong solution for (jl.lj) - (jl.3p goes infinity (for fixed 

initial density) when uo (resp. do) tends to in D A J ,V (resp. B qjP * ). We now aim 
at stating that the system has indeed a global strong solution if uo and do are suitably 
small. This will give Theorem 12.21 

Let (p, u, P, d) be the strong solution given by Theorem 12.11 For any £ > 0, define 

G ,2, C (t) := (||Vpouo|| L 2 + ||Vdo|b) (l + (2(t)h^ and G , 2 := G , 2)0 (t). 
By Lemma 15-H for t < T*, we have 

||^u|| i?(i2) <CGf ^(l + (^t)5) (5.20) 
1 p 

and 

||(Vpu)(t)|| L2 + ||Vd(t)|| L2 < Ce-^G 02 ±(t) 
Hence, starting from ()5.8j) . using (|5.9p . (|5.20j) and the fact that 

ll u • Vu IIl£(L<?) < ll u lli?°(M)l|Vu|| I/ P( i <x.) 

~ C||u|l L oo (D ^l^ l|u|li ?(^ 9 ) ^ C^W' 

||V • ((Vd) T Vd)|| L?(L9) < C||Vd|| ir(i9) ||Ad|| i?(LOO) 

< C|l d ll sfi-i) ||d|| L p (W 3, 9 ) < CG 2 (t), 



u- Vd+ |Vd| 2 d|| L P (L9) 



<C(||, 
we end up with 



<C(||u|| i-i P J|d|| L p (l y2, 9 ) + ||d|| a(i-i) Jl d IU?(w 2 .9)) < CG2 (t), 



G(t) < C((l + £ 2 (t)) + G\t)) + C°(t)G 0>2 tr (1 + (-rit)" 

Once again, the bounds for B p (t) and C p (t) will follow from ()5. 17[) and (j5. 18|) . However, 
in contrast with the previous section, we are going to take advantage of Lemma 15.11 to 
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avoid the appearance of the factor t^~p~^. Indeed, since (L 2 (Q), W 2,q (£l)) td M- W 1,co (£l) 
with i? = y%g, then it follows from Holder's inequality that 

f ||Vu(r)||^dr <C f ||u(r)||J^||u(r)||^„dr 
Jo Jo 

< C f {e-T^G Q2 ±{T)) l -»\\u{T)\\%., q dT 
Jo ' ' P 

<CG l ~ 2 \{t)G\t). 

Now, bounding B p (t) and C p (t) may be done by mimicking the proof of Proposition 15.11 
and we eventually conclude that 

G(t)<Ce °' 2 '^ (l + \\po\\ w i,r) S3 (u (l + tHl + t^)(l + G(t))^)+G 2 (t)j 

(5.21) 

for some positive exponents 63 and £4 depending only on p, q, r. 
Fix a positive T and assume that 

G(t) <8C(l + \\po\\ w i,r) s *U°, te[0,f}. (5.22) 

If the data are so small that 

CGj-^(f)(8C(l + || /90 || W1 ,0 53 ^ ) ,? < In 2, 

then (|5.2ip implies 

G(t) < 2C(1+ llpollwi.O^^Cl + ^U + t^U + GW) 54 ) + G 2 (i)). 
Now, if in addition 

64C 2 (l + ||po||iyi>0 253 ^ < \ and T? (I + T^)^ + 8C(1 + \\p \\ w i, r ) S3 U ) 54 < |, 

then (|5.22p is satisfied with the constant 6C instead of 8C. A standard bootstrap argument 
enables to conclude to the second part of Theorem 12.11 



6. Weak-Strong Uniqueness 

The purpose of this section is to show Weak-Strong Uniqueness in Theorem 12.21 To 
this end, we need to obtain first an energy estimate for the strong solution to system 

flUD-flES). 

Lemma 6.1. Let p,q,r satisfy the same conditions as Theorem \2.1\ and (p, u, P, d) £ 
Mj,q' t be the unique solution to (|l.ip - (|1.3p on Q x [0,To]. Then for any < t < Tq, we 
have 

\ I (p(t)\u(t)\ 2 + \Vd(t)\ 2 )dx + f I (|Vu| 2 + |Ad + |Vd| 2 d| 2 ) dxdr 

^ Jn Jo Jn 2.) 

= - / (p |u | 2 + |Vd | 2 ) dx. 
1 Jn 
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Proof. Integrating (|5.5p over the time interval [0, t], we obtain the energy equality (|6.ip . 

□ 



Now, we proceed to prove weak-strong uniqueness. Let (p, u, II, d) be a global (in time) 
weak solution. On one hand, as the density p satisfies 

(dtp + u- Vp = 
\p|t=o = Po G W^ifl) 

with u G (Lg, c (R + ;fl3(n))) 3 , Theorem 1 in [H] insures that p G C(R+; W 1 ' r ~ (fi)) for all 
r~ < r. On the other hand, we remark that, in view of the regularity of the strong solution 
(p, u, P, d), we deduce from the weak formulation that 

/ pu ■ u dx + / / Vu : Vu dxdr 
Jn Jo Jn 

= / po| u o| 2 ^ x + / / ' (d T u + u • Vu) dxdr — / / (Vd) T Ad • u dxdr 
Jn Jo Jn Jo Jn 



and 



Vd : Vd dx - / |Vd | 2 dx 
Jn 

d • Ad T + u • Vd • Ad - Ad • Ad - |Vd| 2 d • Ad) dxdr 



(6.3) 



o Jn 



Vd| 2 \ ^ JNT/VJ (6.4) 



for almost all t G (0,T ). 
If we write 

pd t M + p\v ■ Vu - Au + VP 

= (p - p)(Stu + u • Vu) + p(u - u) • Vu - V l^f^j ~ (Vd) 1 Ad 

then multiply (|6.4p by u and integrate over f2 x (0, t) to find 

/ / (pd T u + pu ■ Vu) • u dxdr + / / Vu : Vu dxdr 
Jo Jn Jo Jn 

= [ [ ({p - p){d T u + u- Vu) -u + p(u-u) • Vu-u- (Vd) T Ad-u) dxdr 
Jo Jn 

and meanwhile, replace d T d by (jl.lcp in (|6.3p to get 

f Vd : Vd dx- [ |Vd | 2 dx. 
Jn Jn 

2Ad • Ad + u • Vd • Ad + u • Vd • Ad - |Vd| 2 d • Ad - |Vd| 2 d • Ad) dxdr. 

(6.6) 



(6.5) 



o Jn 
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Combining (pHj) . (j6T5|) and (jollj) . we get for almost all t G (0, To), 



/ (pu • u + Vd : Vd) dx + 2 / / (Vu : Vu + Ad • Ad) dxdr 
Jn Jo Jn 

= [ (p |u | 2 + |Vd | 2 ) dx- [ [ (Vd) T Ad • u dxdr 
Jn Jo Jn 



(6.7) 



+ f [ ((p- p){d T u + u • Vu) • u + p(u - u) • Vu ■ u - (Vd) T Ad • u) dxdr 
Jo Jn 

+ [ [ (u • Vd • Ad + u • Vd • Ad - |Vd| 2 d • Ad - |Vd| 2 d • Ad) dxdr. 
Jo Jn 

From (|6.4|) and (|l.lc|) . using the same argument as to get the energy estimate (|6.ip . we 
get for almost all t G (0, To), 

1 [ (p\u\ 2 + |Vd| 2 ) dx + / I (|Vu| 2 + | Ad + |Vd| 2 d| 2 ) dxdr 

2 Jn Jo Jn 

= - f (po|u | 2 + |Vd | 2 )(ix+ / / ((p-p)(a T u + u- Vu) -u + p(u-u) • Vu-u)dxdr. 
^ Jn Jo Jn 

(6.8) 

Here we have used the fact that dtp = — u • Vp. 

Then, adding f)2.1|) and (|6.8p and substracting (j6.7j) . together with the fact that (|5.6p 
and (p, u, n, d) is a weak solution, we obtain 

\ f (p\u - u| 2 + | Vd - Vd| 2 ) dx + f [ (|Vu - Vu| 2 + | Ad - Ad| 2 ) dxdr 
2 Jn Jo Jn 



— ((P ~ / 3 )(^ u + u ' • (u — u) — p(u — u) • Vu • (u — u)) dxdr 



o Jn 

t 



o Jn 



((Vd - Vd) • Ad • (u - u) - u • (Vd - Vd) • (Ad - Ad) 
+ (|Vd| 2 d - |Vd| 2 d) • (Ad - Ad)) dxdr. 
Hence, for almost all t £ (0, To) and for all e > 0, we have 

- f (p\u - u| 2 + | Vd - Vd| 2 ) dx + f [ (|Vu - Vii| 2 + | Ad - Ad| 2 ) dxdr 
2 Jn Jo Jn 

< [ (C e \\d T u + u- Vu||| 3 ||p-p||2 2 +e || u -ii|| 2 6 ) dr + e [ [ \ Ad - Ad| 2 dxdr 
Jo Jo Jn 

+ C £ [ (||Vd||£oo f |d-d| 2 dx) dr+ [ (||Vu|| L oc f p\u-u\ 2 dx) dr 
Jo Jn Jo Jn 

+ ^C e (||Vd + Vd|||oc + ||u|||oo + ||Ad|| 2 3 ) J |Vd - Vd| 2 dxj dr. 

(6.9) 

Here we have used Holder's inequality and Cauchy's inequality with e. 
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Now, we wish to estimate \\p — P\\lz and ||d — d\\ L 2 . We write 

d t (p - p) + V(p - p) ■ u = (u - u) • Vp, (6.10) 

and 

d t {d - d) + u • V(d - d) + (u - u) • Vd = Ad - Ad + |Vd| 2 d - |Vd| 2 d. (6.11) 
Multiply (I6J01) (resp. (IBTTTD ) by p — p (resp. d — d) and integrate over Q x (0, t), we have 

~ [ \p- p\ 2 dx= [ I (p - p)(u - u) • Vp dxdr, 
z Jn Jo Jn 



and 

1 

2 



|d - d| 2 dx 
r-t r ,-•/ 



= -/ /(u-u) ■ Vd- (d-d) dxdr- / [ \Vd - Vd\ 2 dxdr 

+ / / |Vd| 2 |d - d| 2 dxdr + / / (Vd + Vd) : (Vd - Vd) d • (d - d) dxdr. 
Jo Jn Jo Jn 

Employing the same argument as (|6.9p . we get 

\ [ \p-p\ 2 dx< AcellVplllallp-plll.+ellu-Cill 2 ,)^, (6.12) 
z J n io 



|d-d| 2 dx 

n 

t / 1 \ r rt 



< ^C £ ||Vd|| 2 3 + ||Vd||ioo +^\ Jjd-d\ 2 dxdr + ej^ ||u-ii|| 2 6 dr ( 6 13 ) 



+ / [1 + l|Vd + Vd|li °° ] / I Vd - Vdl 2 ,/x,/r. 



Using Sobolev's inequality ||u — u||i6 < C||Vu — Vu|| i 2, we eventually get, from (|6.9 
(16J2D and (I6d3l) . for almost all t G (0,T ), 



/ 

Jn 



u - Ci| 2 + |Vd - Vd| 2 + |p - p| 2 + |d - d| 2 ) dx 

n 

+ [ [ (|Vu - Vu| 2 + | Ad - Ad| 2 ) dxdr 
Jo Jn 



I [ (c £ (t){\p - p| 2 + |vd - vd| 2 + |d - d| 2 ) 

Jo Jn v 

+ C(t)(p\u-u\ 2 + |Vd- Vd| 2 + |d-d| 2 )) dxdr, 



< 



where C e (-),C(-) denote various non-negative measurable functions in L 1 (0, To) occurred 
when we applied the parabolic type estimates for quasi-linear equations (cf. [22], Chapter 
VI, Section 2) to (jl.lcp to obtain d(-,i) 6 C 1,a for some a > 0, with the C 1,a norm 
independent of t. We hence conclude that u = u, d = d and p = p a.e. in x (0,To), by 
applying Gronwall's inequality. 

The proof of Theorem 12.21 is complete. 
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